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We present an atomistic study of the strain field, the one-particle electronic spectrum and the 
oscillator strength of the fundamental optical transition in chemically disordered In^Gai-^As pyra- 
midal quantum dots (QDs). Interdiffusion across the interfaces of an originally "pure" InAs dot 
buried in a GaAs matrix is simulated through a simple model, leading to atomic configurations 
where the abrupt hetero- interfaces are replaced by a spatially inhomogeneous composition profile x. 
Structural relaxation and the strain field calculations are performed through the Keating valence 
force field (VFF) model, while the electronic and optical properties are determined within the em- 
pirical tight-binding (ETB) approach. We analyze the relative impact of two different aspects of 
the chemical disorder, namely: (i) the effect of the strain relief inside the QD, and (ii) the purely 
chemical effect due to the group-Ill atomic species interdiffusion. We find that these effects may 
be quantitatively comparable, significantly affecting the electronic and optical properties of the dot. 
Our results are discussed in comparison with recent luminescence studies of intermixed QDs. 
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I. INTRODUCTION 

Nanometer-size semiconductor quantum dots (QDs) 
have been the subject of many studies in the past years, 
due to their potential applications in optoelectronic de- 
vices and to their peculiar physical properties^ As one 
particularly attractive feature, they offer the possibility 
to tailor the character of the QD electron (or hole) en- 
ergy levels and of the energy of the fundamental opti- 
cal transition by controlling the size, shape and com- 
position of the QD through the growth process. Ex- 
perimentally, InAs QDs in GaAs have been grown both 
by molecular beam epitaxy and metal-organic chemical 
vapor deposition. In most growth processes, nonuni- 
form Ga incorporation in nominally InAs QDs has been 
reported, 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 . 11 ! 12 ! 13 ! 14 ! 15 ! 16 Photoluminescence 
studies of annealed QDs have shown a blue-shift of their 
emission line^iSiiLi 8 ^^ which was suggested to reflect dif- 
fusion of Ga atoms from the matrix material into the 
QD during annealing. However, it is not clear to which 
extent the blue-shift is a consequence of chemical sub- 
stitution (bulk GaAs has a wider band gap than InAs), 
and to which extent it is due to reduced strain in the 
QD after Ga interdiffusion, which also causes a band-gap 
widening. The recently observed change in the photo- 
luminescence polarization anisotropy upon annealing 20 
represents a further interesting but not yet fully under- 
stood QD property. 

From a theoretical point of view, a realistic treatment 
of elastic, electronic and optical properties of such het- 
erostructures must consider a non-uniform In^Gai^^As 
composition profile inside the QD, which we refer here 



as chemical disorder. Several theoretical works deal with 
chemical disorder, either from a macroscopic continuum 
approach, or within a microscopic model. Microscopic 
models provide an atomistic treatment, as required for a 
more reliable description of disordered heterostructures, 
taking into account the underlying zinc-blende structure 
and thus the correct C^v symmetry of pyramidal QDs. 21 
For the elastic properties, previously adopted macro- 
scopic approaches involve a finite element analysis^ 2 ^ or 
a Green's function method^ both in the framework of 
the continuum elasticity theory. Microscopic approaches 
rely on empiric interatomic potentials, such as the Ter- 
soff type^ adopted for truncated pyramidal QDs^ and 
the KeatingSZiS 8 . valence force field (VFF) model, used in 
the study of truncated conical QDsS 

A physical aspect indissociable from atomic interdiffu- 
sion is the strain relief mechanism due to the presence 
of chemical disorder, an effect that has not been high- 
lighted by the previous theoretical studies. We study 
here square-based pyramidal In^Gai-^As QDs within a 
combination of VFF and empirical tight-binding (ETB) 
models, where we distinguish between two different as- 
pects of the chemical disorder on the electronic and opti- 
cal properties, namely the effect of the strain relief inside 
and around the QD and the purely chemical effect due 
to the presence of new atomic species (Ga atoms) pene- 
trating inside the QD. 

From the structural point of view, we calculate the 
strain field inside and around the dot and directly com- 
pare these results with those from a pure InAs/GaAs QD 
of the same size and geometry. This allows a quantitative 
analysis of the strain relief mechanism due to alloying. To 
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simulate the chemical disorder, we employ an atomistic 
diffusion model, where the degree of interdiffusion (and 
thus the degree of chemical disorder) can be controlled, 
so that a direct comparison between a chemically pure 
InAs/GaAs QD and chemically disordered In^Gai-zAs 
dots can be made. 

Regarding the electronic properties, previous studies 
relied on macroscopic approaches such as the single band 
effective mass approximationi2i22i2ii22 or the multiband 
k • p model or on microscopic approaches as 

the empirical pseudopotential model22i2L2£i22i or the em- 
pirical tight-binding (ETB) model. ^ The macroscopic 
models, working with envelope wavefunctions, are ap- 
plicable to smooth composition gradings only^ii and 
cannot properly address the effect of microscopic compo- 
sition fluctuations, which are characteristic of annealed 
samples. We show here that, within ETB, it is possible to 
examine separately how two different aspects of chemical 
disorder affect the QD electronic and optical properties, 
namely the effect of the strain relief inside the QD and 
the purely chemical effect due to In <-> Ga interdiffusion. 
We decouple these effects by performing two indepen- 
dent calculations of the single particle electronic bound 
states and the fundamental optical transition: One in a 
"physical" (strained) QD, and the other in an artificially 
strain- unaffected QD, where only chemical disorder ef- 
fects play a role. Piezoelectric effects were not included 
here, since they become important only for larger QDs. 42 

This paper is organized as follows: In Sec. II we 
present the diffusion model employed to simulate the 
chemical disorder, and we outline the procedure for the 
calculation of the electronic and optical properties within 
the ETB model. In Sec. Ill we present our results, and 
in Sec. IV a summary and conclusions. 



II. FORMALISM 
A. Structural properties 

We start with a square-based pyramidal InAs QD 
with {101} facets and a one-monolayer thick InAs wet- 
ting layer, all embedded in a GaAs matrix. We re- 
strict ourselves for the present purpose to this simple QD 
shape since the relation between the blue-shift and de- 
gree of interdiffusion was found to be only weakly shape- 
dependent^i The pyramid base is 6.8 nm, the height is 
3.4 nm, and the external dimensions of the GaAs matrix 
are 25a x 25a x 17.067a, where a = 5.653 A is the lat- 
tice constant of bulk GaAs. The system contains 85000 
atoms, and periodic boundary conditions are applied. 
Chemical disorder is introduced in the system by allow- 
ing the interdiffusion of Ga and In atoms across the QD 
boundaries. Since the anion sublattice occupation is not 
affected by disorder, we discuss the model in terms of the 
group-Ill species fee sublattice. Interdiffusion is modeled 
atomistically, i.e., each In atom may exchange its position 
with one of its Ga nearest neighbors according to a prob- 



ability p proportional to the concentration of Ga atoms 
around the In atom (p = iVca/12, where Nq^ is the num- 
ber of Ga atoms among its 12 fee nearest neighbors). If 
an exchange takes place, the affected Ga atom is picked 
randomly among the Ga nearest neighbors. We stress 
that the microscopic rules employed to model diffusion 
are compatible with Fick's law of chemical diffusion on 
the macroscopic level. In our diffusion model, one era of 
duration At is completed after all cations in the system 
have been attempted to move. The interdiffusion pro- 
cess is iterated for a discrete number t of eras, and the 
resulting final atomic configuration at t = rAt defines 
the QD to be analyzed. The parameter r quantifies the 
extent of alloying in the system, and simulates the an- 
neal temperature in controlled intermixing experiments. 6 
In order to give some insight about the overall behavior 
to be expected from our assumptions, we present ini- 
tially a description for the evolution of the average oc- 
cupation probabilities at each site. More explicitly, we 
call Pi n (Rj , i) the probability of having an In atom in 
a cation lattice site at the position R^ at a given time 
step t (t = 0, 1, 2, 3, . . . , t). This probability defines the 
average local concentration x of In atoms. Obviously, 
the probability of having a Ga atom at the same position 
and at the same time step is PQ&(R,i, t) = 1 — P[ n (Ri, t). 
The average spatial and temporal evolution of P\ n (Rj , t) 
is described by the equation 

Pi^Ri,*) = Pi„(R l ,t-l) (1) 
1 12 

+ — Pca(Ri, t Pl "( R * + 1 - l ) 

1 12 

- — P In (R i ,t-l)-^P Ga (R i + ^,t-l), 
i=i 

where £j is the j-th nearest neighbor position-vector in 
the fee sublattice. The following points should be men- 
tioned: 

1. In and Ga atoms are treated symmetrically, thus 
the evolution of Pc a (Rj, t) is given by an equation 
analogous to |Q , where the roles of In and Ga are 
interchanged. It follows that the diffusion of In 
atoms into a GaAs-rich region proceeds identically 
to the diffusion of Ga atoms into an InAs-rich re- 
gion. 

2. Ga (In) atoms can penetrate at most r lattice con- 
stants into the QD (into the matrix), and r = 
corresponds to no interdiffusion taking place. 

3. The global concentration does not vary, i.e., the 
total number of cations of each species (In or Ga) 
in the system remains constant. 

A VFF model, parameterized as described in Refs.l2ll 
l43l is then applied to determine the atomic relaxations 
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that minimize the total elastic energy for the given dis- 
tribution of species. In the minimization process, each 
atom is moved along the direction of the force acting on 
it, and the procedure is iterated until the force in each 
atom is smaller than 10 _3 eV/A. 

B. Electronic and optical properties 

The electronic and optical properties are studied 
within an ETB method, adopting a sp 3 s* parametriza- 
tion with interactions up to second nearest neighbors 
and spin-orbit coupling44 Strain effects are included by 
considering both bond length and bond angle deviations 
from ideal bulk InAs and GaAs>^2, Bond length deviations 
with respect to the bulk equilibrium distances affect 
the ETB Hamiltonian off-diagonal elements Vki as 



x (aVR'|i?|acrR) (R'-R)-e , (5) 

where |acrR) represents a general ETB basis vector (a 
runs over the s,p x ,p y ,p z and s* type ETB orbitals, a la- 
bels the spins, R the atomic sites), and C^J R and C^, R , 
are the expansion coefficients of the hole and electron QD 
bound states in the ETB basis. 

Expression (JjjJ can be easily evaluated, since it involves 
all known quantities. Similarly to the electronic prop- 
erties, for the optical properties the strain effects may 
also be completely removed from the calculation. This 
is easily done by using in Eq. JSJ the strain-unaffected 
ETB Hamiltonian matrix elements for (a'er'R'|if |a<rR), 
the strain-unaffected wave function expansion coefficients 
for Cf^L and , RI , and the ideal (bulk) zinc-blende 
interatomic vectors (R — R'). 



Vu (|Ri — Rj 
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where |Rj — Rj| is the actual bond- length and Vki(d^) 
is the bulk matrix element as given in Ref. |44j (k and I 
label the different matrix elements). The exponent n is a 
parameter determined to reproduce the volume deforma- 
tion potentials of InAs and GaAs, whose value was previ- 
ously determined^ as n — 3.40 for all k and I. Strain ef- 
fects may be easily removed from the ETB Hamiltonian. 
The effects of the bond length deformations arc com- 
pletely removed from the Hamiltonian by taking n = 
in Eq. Q . An equivalent transformation causes the effect 
of bond angle deviations from the ideal tetrahedral an- 
gles to be eliminated from the ETB Hamiltonian. Single 
bound hole states \h) and electron bound bound states 
|e) are calculated as eigenvectors of the ETB Hamilto- 
nian, using the folded spectrum method4^ The optical 
transitions in the QD, treated within the electric dipole 
approximation, are quantified in terms of the dimension- 
less oscillator strength 
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where \h) is the initial QD hole bound state, |e) is the 
final QD electron bound state, fuv e h is the transition en- 
ergy, m is the free electron mass, and e is the polarization 
unit- vector. Within ETB, the electron and hole states are 
given by 
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and the electric dipole transition matrix element (e|p ■ 
e\h) can be approximately written as^£ 



(e\p-e\h) ^ 
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III. RESULTS 



Strain field 



Fig. n shows a comparison of the average strain field 
and the local In average concentration between the chem- 
ically pure QD (corresponding to r = 0), given by the 
dotted lines, and the chemically disordered QD (chosen 
here with r = 6), given by the solid lines. For the size of 
QD in this study, r = 6 allows for all but the innermost 
In atoms of the QD to diffuse out. For disordered QD's, 
the results in Fig. ^for each property were obtained by 
averaging over those calculated for an ensemble of 10 dif- 
ferent simulation supercells, all corresponding to r = 6, 
but generated from different sequences of random num- 
bers at each interdiffusion step. In this way, the effect of 
composition fluctuations around the average values given 
in Eq. is reduced. 

The panels on the left side show the xx component 
[panel (a)], the zz component [panel (b)] and the trace 
[panel (c)] of the local strain tensor, as well as the con- 
centration x of In atoms [panel (d)] , along a line oriented 
along the [001] direction and passing through the tip of 
the pyramid. The panels on the right side [(e) - (h)] show 
the corresponding quantities calculated along a line in the 
[110] direction and intersecting the [001]-oriented pyra- 
mid axis at height /i/3 from the base of the pyramid, 
where h is the pyramid height. We observe from frames 
(d) and (h) that, according to our interdiffusion model, 
t = 6 corresponds to a penetration of the Ga atoms in- 
side the QD (and consequently of the In atoms inside 
the GaAs matrix) of about 6 monolayers, i.e. about 17 
A. The error bars shown in the figure indicate standard 

deviations Ae,; 



7<4->-< c tf>7(W-l), where N=W. 



tR 



From the figure we may conclude that 

1. Chemical disorder significantly reduces the abso- 
lute value of the strain field in the regions directly 
affected by the diffusion process, in agreement with 
experimental resultsii On the other hand, very 
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FIG. 1: (Color online) Comparison between the components of the local strain tensor [panels (a)-(c) and (e)-(g)] and the 
concentration x of In atoms [panels (d) and (h)] in a QD of pure InAs (Pure) and in a chemically disordered QD (Disordered) 
with t = 6. The panels on the left side show results calculated along a line oriented along the [001] direction and passing 
through the tip of the pyramid (the value in the horizontal axis corresponds to the position of the wetting layer). The 
panels on the right side show results calculated along a line oriented along the [110] direction and intersecting the [001]-oriented 
pyramid axis at height ft/3 from the base of the pyramid, where h is the pyramid height. The error bars indicate standard 
deviations Ae^ = y/ - (^j) 2 /(N - 1), where AT=10. 
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small changes in the strain field occur in the re- 
gions not affected by interdiffusion, i.e. in the core 
of the pyramid and in the GaAs matrix, at large 
distances from the dot. 

2. If interdiffusion takes place, the strain field varies 
more smoothly than in the case of a chemically pure 
QD. This is a direct consequence of the smooth 
variation of the concentration of In atoms across 
the heterointerfaces of the disordered dots. 



B. Electronic and optical properties 

Fig. [21 shows the calculated eigenenergies of the QD 
bound states as a function of the degree of chemical dis- 
order (characterized by the parameter r). The first two 
electron states ( | el) and | e2) ) are represented in the up- 
per panel, while the first two hole states are shown in 
the lower panel. A chemically pure QD corresponds to 
t = 0. The dashed horizontal lines represent the en- 
ergies of the GaAs bulk conduction (upper panel) and 
valence (lower panel) band edges, delimiting approxi- 
mately the energy range where a QD state is bound. 
The figure shows that the electron state energies increase 
with increasing chemical disorder, while the hole state 
energies decrease, in agreement with previous empiri- 
cal pseudopotential calculations 2 ^. This behavior results 
in an increase of the frequency of the optical emission 
(blueshift), a phenomenon which has been experimen- 
tally observediiZii2ii2i2S The figure shows that, for r = 6, 
the QD gap is about 7% larger than for r = 0. 

Chemical disorder contributes to the results of Fig. [21 
in two ways, namely by the strain relief around the QD 
interfaces (see Fig. QJ, and by the chemical effect due 
to the presence of Ga atoms inside the QD. These two 
effects can be decoupled by comparing the bound state 
energies of a strained QD with those of an artificially 
strain- unaffected QD, as a function of the degree of dis- 
order. Fig. [21 shows such comparison for the energy of 
the electron ground state |el) (upper panel) and of the 
hole ground state \hl) (lower panel). On each panel, the 
uppermost dashed line is a guide for the eye, parallel to 
the QD strain-unaffected energy curve and starting from 
the r = result for the physical QD. The strain relief 
contribution (represented by the solid arrow) can be di- 
rectly compared with the purely chemical effect of the 
disorder, represented by the dashed arrow. We see that 
these two effects are comparable, contributing in oppo- 
site directions for the electron state, and in the same 
direction for the hole state. The purely chemical effect 
can be easily understood: As the interdiffusion increases, 
the concentration x of In atoms in the inhomogeneous 
alloy In-rGai-^As inside the QD decreases. The increase 
(decrease) of the electron (hole) bound state energy as x 
decreases is an alloying effect, so that the electron (hole) 
state energy tends (for x — > 0) to the bulk GaAs conduc- 
tion band minimum (valence band maximum). Results 
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FIG. 2: First two QD bound electron (upper panel) and 
hole (lower panel) state energies as a function of the degree 
of chemical disorder. The dashed horizontal lines represent 
the energies of the GaAs bulk conduction (upper panel) and 
valence (lower panel) band edges, delimiting approximately 
the energy range where a QD state is bound. 



in Fig. [21 show that the chemical effects of disorder are 
partially canceled (enhanced) by the strain relief contri- 
bution for the electron (hole) state. 

We now address the optical properties, focusing on the 
fundamental transition \h\) — ► |el). In Table [I] we com- 
pare the results for a chemically pure QD (interdiffusion 
= off) with those for a chemically disordered QD (in- 
terdiffusion = on) with t — 6. For both cases, an addi- 
tional comparison is made between a strained QD (strain 
= on) and an artificially strain-unaffected QD (strain = 
off). On the first two lines, we show the charge fraction 

AQ = / | ip(r) | 2 d 3 r inside the QD for both the ground 
Jqd 

electron state |el) and the hole ground state \hl). For 
the calculation of AQ in the chemically disordered case, 
the QD border is taken the same as in the chemically 
pure case. The third line shows the oscillator strength 
Jqd of the transition \h\) — ► | el) in the QD for unpolar- 
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TABLE I: Comparison of optical properties between different 
strain states ["strain on" (= "physical" QD) or "strain off" 
(= artificially strain-unaffected QD)] and different degrees of 
interdiffusion ("interdiffusion on" (= chemically disordered 
QD) or "interdiffusion off" (chemically pure QD)). The first 
two lines show the charge fractions within the QD correspond- 
ing to the ground electron (AQ| el )) and hole (AQ| hl )) state. 
The third line shows the oscillator strength Jqd of the fun- 
damental transition \hl) — > |el) in the QD, normalized with 
respect to the oscillator strength /i n As of the fundamental 
transition in bulk InAs. The fourth line gives the degree of 
anisotropy I (Eq. @) of the fundamental optical transition 
with respect to the light polarization direction within the QD 
basal plane. The fifth line shows the oscillator strength /[ooi] 
of the QD fundamental transition for light linearly polarized 
along the [001] direction, normalized with respect to /qd- 
The last line gives the relative change (g — go) /go of the op- 
tical QD gap g with respect to the gap go corresponding to 
the case "strain off" and "interdiffusion off". In the chemi- 
cally disordered case, the error bar (when not negligible) was 
obtained by the same statistical analysis adopted in Fig. Q 



FIG. 3: QD ground electron (upper panel) and hole (lower 
panel) state energy as a function of the degree of chemical dis- 
order. In each panel, we compare the results for the "physical" 
strained QD (QD) with those corresponding to the artificially 
strain-unaffected QD (Strain- unaffected). In each panel, the 
uppermost dashed line is a guide for the eye parallel to the 
QD strain-unaffected energy curve, so that the strain relief 
contribution (represented by the solid arrow) can be visual- 
ized and directly compared with the purely chemical effect of 
the disorder, represented by the dashed arrow. 



ized light, normalized to the oscillator strength /i n As of 
the fundamental transition in bulk InAs for unpolarized 
light. The fourth line gives the degree of anisotropy / 
of the QD fundamental transition with respect to light 
polarization within the pyramid basal plane, denned as 

|( e l|p.e + |M)| 2 -|(el|p-e_|M)| 2 
|(el|p-e + |W)| 2 + |(el|p.e_|/il)| 2 ' { ' 

where e + and e_ are unitary vectors along the inequiv- 
alent basal plane directions [110] and [110] respectively. 
The fifth line shows the oscillator strength /[ooi] of the 
QD fundamental transition for light linearly polarized 
along the [001] direction, normalized with respect to /qd- 
Finally, the last line of Table [I] gives the relative change 



{g—go)/go of the optical QD gap g with respect to the gap 
go corresponding to the case "strain off" and "interdiffu- 
sion off" . As for the case of the electronic properties, the 
direct comparison of the "physical" results with those of 
the disordered case and of the strain-unaffected case al- 
lows us to distinguish between the strain relief effect and 
the chemical effect, both of which are due to chemical 
disorder. 

From the results in Table [IJ we arrive at the following 
conclusions: 

1. The first two lines show that chemical disorder 
reduces the confinement of the QD bound states 
through the partial relief of the strain field, while 
the chemical effect does not directly contribute. In 
fact, chemical disorder reduces the charge fractions 
AQ| e i) and AQ\hi), while no changes are observed 
for the strain-unaffected calculation. The smaller 
confinement of the QD bound state wave functions 
in the chemically disordered case is consistent with 
the results of Fig. [3 where all electron and hole 
bound states become shallower when chemical dis- 
order increases. 

2. The third line indicates that chemical disorder sig- 
nificantly enhances (by about 40%, in the case con- 
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sidered here) the oscillator strength /q£> of the fun- 
damental optical transition, in qualitative agree- 
ment with experimental resultsii2*i& This effect is 
primarily due to the modification of the strain field 
due to chemical disorder, because in the strain- 
unaffected case Jqd does not significantly vary. 

3. From the fourth line, we observe that chemical dis- 
order strongly reduces the in-plane asymmetry J, 
in accordance with previous experimental results^ 
This is a direct consequence of the partial relief of 
strain due to disorder. In fact ^ in a pyramidal QD 
the asymmetry of the oscillator strength of the fun- 
damental optical transition between the directions 
[110] and [110] is a direct consequence of asym- 
metry of the strain field between these directions, 
which is in turn a consequence of the Ci v symme- 
try. This can be deduced observing that in the 
strain-unaffected case / vanishes. This result could 
be experimentally exploited to detect, among dif- 
ferent samples containing QDs of similar geometry, 
those characterized by the higher chemical purity. 
In fact, these samples will be those with the higher 
asymmetry of the absorption coefficient of the fun- 
damental optical transition (which is proportional 
to /), for in-plane polarized light. 

4. The fifth line of Table [i] implies that chemical dis- 
order weakens the fundamental optical transition 
for perpendicularly polarized light. This is a con- 
sequence of the strain relief inside the QD: In the 
limit of complete relief (QD strain-unaffected) this 
transition is strictly forbidden4& 

5. The last line summarizes the different effects con- 
tributing to the blue shift in the fundamental opti- 
cal transition with respect to a hypothetical tran- 
sition energy go where both effects are removed. 
We see that strain and chemical disorder increase 
the QD gap by the same order of magnitude. We 
note that calculations for the relative blue shift pre- 
sented in Ref . systematically underestimate this 
quantity as compared to the experimental results in 
Ref. (see Fig. 7 in Ref. l4l|) . This discrepancy is 
probably due to the simplified theoretical descrip- 
tion adopted there, where strain effects were not 
taken into account. 

Finally, we analyzed the z-component of the built-in 
dipole moment of the electron-hole pair, and how disor- 
der affects it. Such dipole moment experimentally shows 
up as a Stark shift of the emitted light from a QD- 
LED under applied electrical fields For pure pyramidal 
InAs/GaAs QDs, this dipole moment points towards the 
base of the pyramid, i.e. the center of mass of the electron 
ground state lies above that of the hole ground stated 
However, in the case of truncated pyramidal In^Gai-ajAs 
QDs, with x increasing from the base to the tip of the 
pyramid, the dipole moment may have an opposite ori- 
entation, i.e. the center of mass of the hole state can 



sit above that of the electron stated Some authors have 
argued that such inversion occurs also for QDs having an 
In-rich core with an inverted-cone shape. This inverted- 
cone shape has been observed in truncated-cone nomi- 
nal Ino.8Gao.2As QDsS and Ino.5Gao.5As QDs^i In our 
case, the dipole moment is always directed towards the 
base of the pyramid, i.e. the electron ground state sits 
always above the hole ground state, both for the pure 
and the disordered QD. This is because we have neither 
a truncated pyramidal shape nor an In-concentration in- 
creasing from the base to the tip of the pyramid (see 
Fig. However, we observe that the disorder decreases 
the dipole moment of the dot. In fact, in the strained dis- 
ordered case, the center of mass of the electron state lies 
2.8 A above that of the hole state, while in the strained 
pure case this separation is 3.5 A. 



IV. SUMMARY AND CONCLUSIONS 

We presented an atomistic interdiffusion model to sim- 
ulate the composition profile of chemically disordered 
pyramidal In^Gai-^As QDs buried in GaAs matrices. 
Calculations for the strain field inside and around the dis- 
ordered QDs were compared to the strain field of chem- 
ically pure InAs QDs, showing that chemical disorder 
significantly reduces the absolute value of the strain field 
inside the QD, giving rise to smoother variations of this 
field across the heterointerfaces. Furthermore, we ana- 
lyzed the consequences of chemical disorder for the elec- 
tronic and optical properties within an ETB model. Our 
treatment allowed us to distinguish between two effects 
of the chemical disorder, namely the relief of the strain 
inside the QD, and the purely chemical effect due to the 
presence of new atomic species (Ga atoms) penetrating 
inside the QD. We showed that these two components 
of disorder have comparable effects on the QD electronic 
spectrum, while for the optical properties the strain relief 
effects are more relevant. In particular, we showed that 
strain relief (i) reduces the charge confinement (inside the 
QD) of the electron and hole bound state wave functions, 
(ii) significantly enhances the oscillator strength of the 
fundamental optical transition, (iii) strongly reduces the 
asymmetry of the oscillator strength of the fundamental 
optical transition between the directions [110] and [110] 
for in- plane polarized light, and (iv) strongly reduces the 
oscillator strength of the fundamental optical transition 
for perpendicularly polarized light. 

Our results help to explain experimental findings for 
the optical properties of intermixed InAs/GaAs QDs. 
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